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ABSTRACT 



Context. The availability of precisely determined frequencies of radial and non-radial oscillation modes in red giants is finally paving 
the way for detailed studies of the internal structure of these stars. 

Aims. We look for the seismic signature of regions of sharp structure variation in the internal structure of the CoRoT target HR7349. 
Methods. We analyse the frequency dependence of the large frequency separation and second frequency differences, as well as the 
behaviour of the large frequency separation obtained with the envelope auto-correlation function. 

Results. We find evidence for a periodic component in the oscillation frequencies, i.e. the seismic signature of a sharp structure 
variation in HR7349. In a comparison with stellar models we interpret this feature as caused by a local depression of the sound speed 
that occurs in the helium second-ionization region. Using solely seismic constraints this allows us to estimate the mass (M = 1.2^ 
Mq) and radius (R = 12.2^-1 R Q ) of HR7349, which agrees with the location of the star in an HR diagram. 
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1. Introduction 

Red giant stars are cool, highly luminous stars. They play 
a predominant role in stellar, galactic, and extragalactic as- 
trophysics because they serve as distance and age indicators 
for g l obular clusters and ne arby galaxies (see e.g. iLee et al.1 
1993; iGirardi & Sa laris 2001). Moreover their observed chem- 
ical composition is a fundamental ingredient in t he study of the 
chemical evolution of their host galaxy (see e.g. ICatelanl [20091 
and references therein). 

All red giants have an extended and diluted convective enve- 
lope surrounding a dense core, which makes their structure and 
therefore their pulsation properties very different from that of 
our Sun. As in solar-like stars however, their convective enve- 
lope can stochastically excite acoustic oscillation modes. While 
stochastic oscillations h ave been firmly detect ed in a few bright 
G and K giants since IFrandsen et al.l d2002l) . a major break- 
through was achieved thanks to the p hotometric space -based ob- 
servations with the CoRolQ satellite (Baglin et a l.l2006l) . leading 
to the unambiguous detection of radial and non-radial oscilla- 
tion modes in thou s ands of red giants dDe Ridder etafl [2 009 ; 
iHekker et all l2009t ICarrier et all l2010t iMosser et al.l I2016T . 
Moreover, the Kepler space-telescope has also started monitor- 
ing red giants for photometric variations and has already de- 
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tected radial and non-radial oscillation mod es in a large number 
of field giants as w ell as cluster members (Bedding et al. 2010a; 
IStello et at]|2010l) . As the precise detection of oscillation modes 
in the Sun fostered helioseismic inferences in the 1990s, we now 
have the chance to start probing the internal structure of these 
cool evolved stars thanks to the recent advances in observational 
stellar seismology. 

Stellar global acoustic oscillation modes show up in the fre- 
quency spectrum as discrete peaks. While the average spacing 
between the frequencies v of modes of the same spherical de- 
gree (large frequency separat ion, Av) represent s an estimate o f 
the mean density of the star (Vandaku rovll 19671 lTassoulH r980). 
periodic deviations from a constant Av are signatures of so- 
called acoustic glitches, i.e. regions of sharp-structure variation 
in the stellar interior. Any localised region of sharp variation 
of the sound speed c induces in the frequencies an oscillatory 
component with a periodicity related to the sound-t ravel time 
measured from that region to the surface of the star dVorontsovl 
1988; lGoughl [l990). Its amplitude depends on the sharpness of 
the glitch and decays with v, because as the frequency increases, 
the wavelength of the mode becomes comparable with or less 
than the extent of the glitch. These sharp variations are found 
in zones of rapidly changing chemical composition, in ioniza- 
tion zones of major chemical elements, or in regions where the 
energy transport switches from radiative to convective. 

In the Sun, the detection and interpretation of a large num- 
ber of acoustic oscillation modes led to the determination of 
the depth of both the solar convective envelope and of the he- 
lium second-ionization zone, and allowed us to estimate the he- 
lium abund ance in the envelope, otherw ise inaccessible by other 
means (see lChristensen-D alsgaard 2002, for a review). 
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The possibility and relevance of detecting sharp 
structure variations in other main-sequence stars, as 
well as first attempts to look for these signatures, have 
been addressed in several papers in the literat ure (see 
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Here we report on the detection of a sharp-structure varia- 
tion in the red giant C0R0T target HR7349 (HD181907), a star 
with a physical structure very different from the Sun, yet show- 
ing global oscillations modes that allow us to probe its internal 
structure. 



2. Evidence for an oscillatory component in the 
frequency spacings 

ICarrier et al.l (|2010) have determined the eigenfrequencies of 
HR 7349. They have shown that the first frequency differences 
Ay„ f = v„/ - v„_i.f are modulated. In order to investigate 
this modulation, we have analysed these first frequency differ- 
ences and compared the m to the function Av(v) derived from 
the method developed by iMosser" & Appou rchauxl (12009). This 
envelope autocorrelation function (EACF) is obtained from the 
auto-correlation of the times series, performed as the Fourier 
spectrum of the windowed Fourier spectrum; the windows are 
Hanning filters covering the frequency range where significant 
excess power is observed. Because the modulation is rapid, we 
had to use narrow filters, with a FWHM varying between 1 and 
2 times the mean value of the large separation. An EACF com- 
puted with a narrow filter benefits from the low filtering, but suf- 
fers from larger error bars. 

Fig. Q] shows the large separations for I = 0, 1, and 2 super- 
imposed on the function Av(v) obtained with a filter width equal 
to 1 .4 times the mean value of the large separation, which rep- 
resents an acceptable compromise between the filtering of the 
signal and the noise level. We remark that both values and error 
bars of Av(v) and Av„j closely agree. The modulation seen in the 
function Av(v) is present in all autocorrelations performed with 
a filter that is narrow enough to respect the Shannon criterion. 

We have detrended the low-frequency gradient of Av(v) in 
order to emphasize the modulation, which is of the order of 
7 /iHz (a period of about To = 0. 14 Ms). 

That the modulation is seen with two independent methods, 
is present over more than 2.5 periods, and remains detectable 
even if a significant level of noise is added to the spectrum makes 
us confident that the oscillation found in the frequency differ- 
ences is genuine. 

The uncertainty of the period To of the oscillatory com- 
ponent was estimated with a Monte Carlo simulation. We as- 
sumed that the frequencies are independent and are Gaussian- 
distributed around their observed value with a standard dev iation 
equal to their l-cr error bar as listed by ICarrier et aj] (l2010h . We 
then generated 5000 realizations of the frequencies, computed 
the frequency differences, and applied a non-linear regression 
procedure. Each time a cost function S , defined as 



S(P) = (y-ii<J3)yv- l (y-nm 
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Fig. 1. Large separations Av„/ computed by [Carrier et alJ {2010) 
(t = 0: black squares; 1: stars; 2: triangles) are superimposed 
on the Av„/ function given by the autocorrelation of the time se- 
ries (solid thick line). The ±l-cr error bars are given by the two 
adjacent curves; they agree with the error bars of Av n j. Both 
Av„/ and Av n j show a periodic modulation with a period of 
the order of 0.14 Ms. The zebra regions indicate 1-Av wide fre- 
quency ranges. The inset shows the filter used for performing 
the envelope autocorrelation function, with a full-width at half- 
maximum 6vu equal to 1.4 times the mean value of the large 
separation. 



differences, and /u(J3) is the matrix containing the non-linear 
model 



fi(j3; v ) = a + bv + A sin(27r T >v) + B cos(2^ To v) 



(2) 



with B = (a, b,A, B, T )) the fit parameters, and v the covariates. 
The uncertainty of To was then derived as the standard devia- 
tion of the sample of To values. As a result, we obtain a period 
To = 0.13 + 0.02 Ms from the component detected in Av and 
To = 0.14 + 0.03 Ms when fitting the behaviour of t he second 
frequency differences A2V (as defined in lGoug h 1990). 

3. Interpretation of the oscillatory component 

As recalled in Sect. [1] the periodicity of the component allows 
us to infer the location of the glitch in terms of its acoustic 
de pth or in terms of its relative acoustic radius (as suggested 
bv lBallot et al1l200l) . For HR7349 we find: t /T = 0.55 + 0.07, 
where t is the sound travel time from the centre to the glitch, and 
T is the acoustic radius of the star (T = (2 (Av)) -1 ). 

To interpret this observational evidence, we consider in Fig. 
|2] the location of acoustic glitches in t he structure of red -giant 
models computed with the code ATON dVentura et al.l 2008). We 
find that at about half the acoustic radius the sound speed has a 
sharp variation owing to a local depression of the first adiabatic 
exponent Ti in the helium second-ionization zone. We there- 
fore interpret the periodic component detected in the frequen- 
cies of HR7349 as the signature of helium second ionization. 
This interpretation is further supported by a direct comparison 
between observed and numerically computed pulsation frequen- 
cies in models of red giants (see Sect. |4j. 

We notice that while in the Sun and solar-like stars the base 
of the convective envelope is located at about half the stellar 
acoustic radius (see Fig. ff]), in red giants it is located at a small 
relative acoustic radius (typically less than 10 _1 for the models 
considered, depending on the mass and evolutionary state). 



was minimized, where y is the matrix of observed frequency dif- 
ferences, V is the covariance matrix of the observed frequency 



2 71 = (d lnp/<91np),, where p, p, and i - denote pressure, density, and 
specific entropy. 
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Fig. 2. a) Behaviour of y\ as a function of t/T in a 1 .2 M , 12 R Q 
model: the region where helium undergoes its second ionization 
induces a local minimum in y\ . In the same model the base of 
the convective envelope is located at t/T = 2 10~ 2 . b) Location 
of acoustic glitches in a model of the Sun. 

Furthermore, we find that in the models the relative acoustic 
radius of the He second-ionization region (tn e n/T) and the av- 
erage Av determine the mass and radius of the star (see Fig. [3). 
We ascribe this dependence to the the fact that the temperature 
stratification in the e nvelope of red-giant stars dep ends mostly 
on M and R (see e.g. lKippenhahn & Weigertlll990l) : a thorough 
discussion on this topic will be reported in detail in a forthcom- 
ing paper. 

On the basis of seismic constraints only, we can thus derive 
from Fig. |3]the mass of HR7349 (M = 12™ M ) and, through 
the scaling of Av with the mean stellar density, its radius (R = 
12.2+^ ^ R Q ). 

Because the parallax of this nearby red giant is known, 
an independent radius estimate of HR7349 can be obtained 
via the precise abs olute luminosity and effective temperature 
(ICarrier et al.l l2010h . leading to R — 12.2 +1.1 R G , in agree- 
ment with the seismically inferred radius. We also notice that 
the mass and radius estimates are also compatible to those ob- 
tained by the combination of the frequency corresp onding to the 
maxim um oscillation power (v max ), (Av), and r e fj (Carrie r et alj 
l20loh . 

4. Direct comparison with numerically computed 
oscillation frequencies 

We consider a red giant model of 1 .2 M Q and 12 R n compute d 
with the stellar evolutionary code ATON (IVentura et al.l l2008). 
The model has an initial helium mass fraction Y = 0.278 and 
a heavier-elements mass fraction Z = 0.02. For this model we 
computed adiabatic oscillation frequencies of spherical degree 
I — 0, 1 and 2 in the domain 15-40 /zHz with an Eulerian version 
of the code LOSC dScuflaire et al.ll2008l) . 

In a red giant model, differently from the Sun, non-radial os- 
cillation modes in the frequency domain of solar-like oscillations 



variation inside a red-giant star 
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Fig. 3. Behaviour of the relative acoustic radius of the second he- 
lium ionization region (fHeii/7") as a function of the average large 
frequency separation in stellar models with different mass (0.7 
to 3 Mq), initial chemical composition (helium mass fraction 
Y = 0.25, 0.278, heavier-elements mass fraction Z = 0.02, 0.01), 
and convection-efficiency parameter (ctmlt = 1.6, 1.9). The po- 
sition of HR7349 in this diagram is represented by the black dot 
with error bars. 

dKieldsen & Beddinil l 19951) are expected to have a mixed pres- 
sure and gravity nature. In the model we considered (see Fig. 
[3} the spectrum of t = 1 and 2 modes is mostly populated by 
modes with a dominant gravity-like behaviour. These modes ex- 
hibit a high inertia compared to radial modes, and therefore show 
surface amplit udes significantly smaller than the purely acoustic 
radial mod es (Dziembowski et al. 2001; Christe nsen-Dalsgaardl 
l2004tlDupret et al.l2009l:lEggenberger et al.l2010l) . Nonetheless, 
among this large number of non-radial oscillation modes, few 
modes are efficiently trapped in the acoustic cavity of the star 
and have inertias similar to those of radial modes (see Fig.|4]i. 

We then compute Av(v) and A2v(v) by selecting radial modes 
and the I = 1 modes corresponding to local minima in the in- 
ertia, hence strongly trapped in the acoustic cavity. We find an 
average value of the large frequency separation of 3.4 yuHz and 
clear periodic components in Av(v) (open squares in Fig. [5j and 
A 2 v(v). 

The comparison between the model and observed frequency 
separations clearly supports the interpretation that the periodic 
component is caused by the local depression in the sound speed 
in the second helium ionization zone. This feature is well repre- 
sented by current models that have a mass and radius determined 
from the periodicity of the component and the average value of 
the large frequency separation (see Fig. [3}. 

We notice that the periodic component related to the He 
second-ionization zone is expected to have an amplitude de- 
caying with frequency, as shown in numerically computed fre- 
quencies (o£eji_sc|UOTes_jn_J^^^ 

mations ( Monte iro & Thompson 1998.llloudek & Goughl2007h . 
Although current error bars are too large to draw a firm state- 
ment, the amplitude of the observed component computed with 
t — and 1 modes seems to show the expected behaviour, while 
we notice that the observed I = 2 modes show a larger scatter. 
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Fig. 4. Mode inertias of radial modes (full dots), I = 1 modes 
(asterisks) and 1-2 modes (open triangles) for a 1.2 M red- 
giant model. 
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Fig. 5. Open squares represent Av„/ computed from I — 0, 1,2 
adiabatic frequencies in a 1.2 M red giant model. The solid line 
shows a sinuso idal component with amplitude decreasing with 
frequency (see Monteiro & Thompson 1998; Houdek & Gough 
2007 ) fitted to th e I — 0, 1 large separation determined by 
Carrie r et al 2010 (dots with error bars). 



5. Conclusions 

In summary, we have shown that the acoustic pulsation modes 
detected by CoRoT in the red giant star HR7349 bear the sig- 
nature of a sharp-structure variation inside the star. Comparison 
with stellar models allows us to interpret this feature as caused 
by a local depression of the sound speed occurring in the helium 
second-ionization region. 

Besides representing the first seismic inference of a local fea- 
ture in the internal structure of an evolved low-mass star, this de- 
tection allows a mass (M = 1.2^ M ) and radius (R = 12.2+fl 
R Q ) estimate based solely on seismic constraints. Moreover, for 
this nearby star, we could also check that our radius estimate is 
compatible with the one based on luminosity and effective tem- 
perature {R = 12.2 +1.1 R Q ). This additional test reinforces the 
proposed interpretation and approach, which could be applied 
to the thousands of pulsating giants of unknown distance that 
are currently being obs erved with the space sate llites CoRoT 
(see Hekker et al.ll2009 | and, in particular , Fig. 6 i n Mosser et al.l 



2010) and Kepler dBorucki et all l2010t iGilliland el "alJ l2010t 
Beddi ng et al. 2010a|). A reliable seismic estimate of the mass 
and radius of these stars wou ld represent an esse ntial ingredient 
for stellar population studies (Mig lio et al.l l2009) and for charac- 
teri zing planets orbiting around these evolved distant stars (see 
e.g. lHatzes & Zechin eister 20Q3). 



We finally recall that the detectability of the signature of He 
ionization can potentially lead us to a seismic estimate of the en- 
velope helium abundance in old stars. Indeed, as shown for the 
Sun and solar-like stars, the amplitude of the periodic component 
depends on the envelope helium abundance (see e.g. Bas u et al.l 
120041: iHoudek & Goughl 120071 and references therein). While 
CoRoT and Kepler observations will provide other targets and 
further reduce the uncertainty on the oscillation frequencies, a 
thorough study on the required precision in terms of seismic and 
non-seismic observational constraints, as well as and in terms of 
models, should be undertaken to aim for reliable seismic esti- 
mate of the envelope helium abundance in giants. 
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